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Abstract 



We study a two-dimensional (2D) spin-half Heisenberg model related to the 
quasi 2D antiferromagnets (Ba, Sr)2Cu$0±Cl2 by means of exact diagonal- 
ization and spin-wave theory. The model consists of two inequivalent in- 
terpenetrating square-lattice Heisenberg antiferromagnets A and B. While 
the antiferromagnetic interaction J aa within the A subsystem is strong the 
coupling Jbb within the B subsystem is much weaker. The coupling Jab be- 
tween A and B subsystems is competing giving rise for interesting frustration 
effects. In dependence of the strength of Jab we find a collinear Neel phase, 
non-collinear states with zero magnetizations as well as canted and collinear 
ferrimagnetic phases with non-zero magnetizations. For not too large values 
of frustration Jab, which correpond to the situation in (Ba, Sr)2Cu2,O^Cl2, 
we have Neel ordering in both subsystems A and B. In the classical limit 
these two Neel states are decoupled. Quantum fluctuations lead to a fluctu- 
ational coupling between both subsystems ('order from disorder') and select 
the collinear structure. For stronger Jab we find evidence for a novel spin 
state with coexisting Neel ordering in the A subsystem and disorder in the B 
subsystem. 



I. INTRODUCTION 



The exciting collective magnetic properties of layered cuprates have attracted much attention 
over the last decade. A lot of activity in this field was stimulated by the possible connection 
of spin fluctuations with the phenomenon of high-temperature superconductivity. But, the 
rather unusual properties of quantum magnets deserve study on their own to gain a deeper 
understanding of these quantum many-body systems. In recent years some of those materials 



such as BCL2CU3O4CI2 and Sr 2 Cu 3 04Cl2 have been studied experimentally and theoretically 
in more detail [1-5]. The most important difference between (Ba, Sr) 2 Cu^04Cl2 and their 
parent compound La^CuO^ is the existence of additional C«(i?)-atoms located at the centre 
of every second Cw(A)-plaquette. Both subsystems A and B form square lattices, however, 
with different orientations and lattice constants. This A-B lattice is illustrated in Fig. la. 
Both copper sites Cu(A) and Cu(B) carry spin half, i.e. quantum fluctuations are important. 
Since the magnetic couplings Jaa between A spins as well as Jbb between B spins are 
antiferromagnetic and the coupling Jab between A and B spins is frustrating we have a 
system of two competing antiferromagnetic spin-half subsystems. 

Three-dimensional examples of two interpenetrating antiferromagnets like garnets 
Mn 3 Cr 2 Ge 3 0i2 or (Fe x Gai- x )2Ca 3 Ge 3 0i 2 were discussed by several authors (see e.g. [6-8]). 
For the quasi two-dimensional cuprates like (Ba, Sr) 2 Cu2,OiCl2 the quantum fluctuations 
are more important than in the three-dimensional garnets and the interplay of competing 
interactions with strong quantum fluctuations may lead to interesting magnetic phenomena. 

Noro et al. [1] reported two magnetic phase transitions at T A = 320K and at T B = 
for Ba 2 Cu 3 4 Cl 2 being attributed to respective antiferromagnetic ordering of the Cu(A) 
and Cu(B) spins. Both critical temperatures differ in one order of magnitude indicating a 
strongly antiferromagnetic coupling between Cu(A) spins and a comparatively small antifer- 
romagnetic coupling between Cu(B) spins, which is confirmed by band-structure calculations 
[9]. According to Chou et al. [3] the weak ferromagnetic moment found experimentally [2] 
could be understood as a consequence of bond-dependent interactions such as pseudodipolar 
couplings. 

The minimal model to describe the main magnetic properties of the competing antiferro- 
magnets on the A-B lattice is the antiferromagnetic Heisenberg model with three exchange 
couplings Jaa, Jbb and Jab- In what follows we call this model A-B model. Some prelimi- 
nary results for a finite system of Af = 24 sites were reported in the conference paper [10]. 
However, to describe the weak ferromagnetism observed in these compounds anisotropic 
interactions seem to be needed [3-5] . 

In this paper we want to study the influence of strong quantum fluctuations and frustration 
on the ground state of the A-B model using spin-wave theory and exact diagonalization. 
The paper is organized as follows: In Section II we introduce the A-B model and illustrate 
the classical magnetic ground-state phases in the considered parameter region. In Section 
III we present an exact-diagonalization study of the ground-state phases and in Section IV 
the linear spin-wave approach is used to analyse the Neel phase realized for small Jab i n 
more detail. In Section V a summary is given. 
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II. THE A-B MODEL AND ITS CLASSICAL GROUND-STATE PHASES 



We consider the Hamiltonian (cf. Fig. lb) 
H = Jaa E 

+ Jbb E + Jab E (1) 

{m£A,neA} (meB,neB) (m£A,n£B) 

where the sums run over neighbouring sites only. Jaa and Jbb denote the antiferromagnetic 
couplings within the A(i?)-subsystems, respectively. We focus our discussion on parameters 
Jaa = 1 and Jbb = 0.1 which corresponds to the situation in (Ba, Sr^CuzO^Clz. The 
value of the frustrating inter-subsystem coupling Jab is less reliably known. We consider 
antiferromagnetic Jab and use it as the free parameter of the model. The lattice consists of 
J\f = 3N spins with three spins per geometrical unit cell and ten couplings in it. 

We start with the discussion of the classical ground state, i.e. the spins S n are considered 
as classical vectors of length s. Varying Jab we have altogether five ground-state phases, 
see table I. Two of them (I and III) have planar spin arrangement, two (II and IV) are 
non-planar and one (V) is collinear. Without loss of generality we choose in this section for 
the description of planar spin ordering the x-y plane. We start from weak inter-subsystem 
coupling J A b ~ 0. Then we have Neel ordering in both subsystems (phase I, table I). 
These two classical Neel states shown in Fig. 2 are decoupled and can rotate freely with 
respect to each other, i.e. the ground state is highly degenerated and this degree of freedom 
is parametrized by the angle ip. The corresponding magnetic unit cell contains six spins. 
(Thus, later in Section IV we have to introduce six different magnons in the spin-wave theory 
for this phase.) 

At Jab = 2V J aaJbb there is a first-order transition from the Neel phase I to the non-planar 
ground-state phase II. Phase II is illustrated in Fig. 3. The corresponding magnetic unit cell 
contains 12 pins and is therefore twice as large as the magnetic unit cell of the Neel phase 
I. In this state we have eight different spin orientations characterized as follows 

S A \ = s i^-^Y cos (o),-^y cos («) » sin (±«)^ » 

s£ = s {^y cos ( a ) » cos ^ ' sin ( ±a ^j ' 
s A [ = s i^-^Y cos ^ ' cos ( a ) ' ~~ sin ( ±a )) ' 
Sa 4 = s (^y cos ( a ) ' ~^y cos ^ ' ~ sin ( ±a )) ' 



Sg =s (-1,0,0), Sg = s (0,-1,0), 

ill _ „ ft n n\ ell 



Sg=s (1,0,0), S^ 4 = s (0,1,0), (2) 
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where a is given by a = arcos [Jab/ Jaa) ■ Obviously neighbouring B spins are per- 
pendicular to each other and consequently the energy Ejj is independent of Jbb (see table 
I). The in-plane xy components of the A spins of neigboring spins are also perpendicular, 
however, there are finite off-plane z components. These off-plane components (proportional 
to sin a, see eq. (2)) decrease with Jab and become zero at Jab — Jaa, i-e. we have a 
second-order transition from the non-planar phase II to the planar phase III at this point. 
The spin orientations of phase III are given by eq. (2), too, but with the additional condition 
a = (see Fig. 3). In phase III neighbouring A-spins as well as neighbouring 5-spins are 
perpendicular to each other and consequently the energy depends on Jab only. All three 
phases I, II, III have a zero net magnetization S tota i = 0. 

Further increasing Jab favours an antiparallel alignment of the A spins relative to the B 
spins and the planar phase III gives way to a non-planar phase IV where the in-plane xy 
components are aligned as in phases II and III (see Fig. 3). The off-plane z components 
in phase IV are given as S* eA = +s cos 6 for A spins and as S* eB = —s y/l — A 2 sin 2 9 for 
B spins. This phase IV can be denoted as canted ferrimagnet [11] and has a net magnetic 
moment S A in the A subsystem and Sb in the B subsystem resulting in a finite total magnetic 
moment S tota i 



S 1 / 



^ = ^^ sN 



1-AW0; SZ al = \S A -S B 



(3) 



with 



sin# = 





1 - 


I Jaa 
\Jab 


1 Jab 
^ 4Jbb 


-v 


24. j aa _ 
Jbb 


t2 72 

h 16^ + 4^ 

J AB J BB 1 


) 


g Jaa 

Jbb 


+ 2 


J 2 

AA _l_ 

~Ti r 

J AB 


*4B_ + 

SJ BB 


f 3Jaa 
\2Jab 


Jab \ 
SJbbJ 1 


/_244 dA + 16 

/ Jbb 


j2 72 
J AA _|_ J AB 

J AB ^BB 



(4) 



and 



Jab 



V2 



Jaa 



3 ^aa +2 ^aa + -Jab 



BB 



T 2 

J AB 



8 Jbb 



2y/2 Jbb Jab 

This total moment increases with Jab- The energy of phase IV is given by 
ElV = \Jaa (l - sin 2 (9 ) + \Jbb (l - A 2 sin 2 ) 



s 2 Af 



Jab + ^(l-sin 2 ^) (l - A 2 sin 2 o fj 



(5) 



(6) 



The phase boundary of the second order phase transition between III and IV is given by 
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J^t IV = V2J AA + Jj ^ + \ ^J\a + 24 J AA J BB + 2 lf B (7) 
and yields J A £~ IV = 3.099 for J AA = 1 and J BB = 0.1. 

Finally, for large J AB the A and B spins are fully polarized along the z axis, i.e. S n< z A = 
(0, 0, +s) and S ne B = (0, 0, — s) and a collinear ferrimagnetic phase V is realized. The phase 
boundary of this second order phase transition between IV and V is given by 

Jab~ V — 2Jaa + Jbb + \J 4J AA + Jbb ( 8 ) 

leading to J 1 /^ = 4.102 for J AA = 1 and J BB = 0.1. 

III. THE QUANTUM GROUND STATE - EXACT DIAGONALIZATION 

To discuss the influence of quantum fluctuations on the classical phases studied in the 
last section we use the Lanczos algorithm to calculate the quantum ground state of the 
Hamiltonian (1) for a finite lattice of M = 24 spins (Fig.4). Again we choose the parameters 
J aa = 1) Jbb = 0.1 appropriate for (Ba, Sr) 2 Cu^O^Cl2 and consider J AB as the free 
parameter. For M = 24 sites we have 16 A spins and 8 B spins. Since the maximal 
magnetic unit cell of the classical ground states contains 12 spins the M = 24 system has 
the full symmetry of the classical ground state in the considered parameter region. To reduce 
the Hilbert space we used all possible translational and point symmetries of the A-B lattice 
as well as spin inversion. The use of the symmetry allows to classify the different quantum 
ground states by their symmetry. 

To compare classical and quantum ground-state phases we present the spin-spin correlations 
in Figs. 5, 6 and 7. The Neel phase I is present also in the quantum model. However, the 
quantum fluctuations lift the classical degeneracy and both subsystems couple. The fluctua- 
tional coupling is known as order from disorder effect [12,7] and selects a collinear quantum 
state with a finite A-B spin correlation in the quantum Neel phase I (see Fig. 7). Moreover, 
the transition to phase II is shifted to higher values of J A b indicating that quantum fluc- 
tuations favour collinear versus non-collinear states (see e.g. [13]). The frustrating coupling 
weakens the A-A and B-B spin correlations in the Neel phase I; this weakening is stronger 
for the B-B correlations than for the A-A correlations. Hence, a disordered quantum ground 
state similar to the J\ — J 2 model [14-18] seems to be possible and will be discussed in more 
detail in next Section. 

The quantum Neel phase gives way to a fairly complex spin state at J A b ~ 1-02 up to 
Jab ~ 3.13. This state is also a singlet S = as the quantum Neel state. The A-A and the 
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A-B spin correlations of the quantum model follow qualitatively the classical curves (see Figs. 
5 and 7). However, we see some jumps in the correlations connected with level crossings of 
ground states belonging to different lattice symmetries. Most likely these level crossings may 
be attributed to finite-size effects. The change of A-A and A-B correlations at Jab ~ 1-02 
is small. However, the B-B correlations change strongly at this point. Contrary to the 
classical model where the nearest-neighbour B-B correlation is zero and the next-nearest- 
neighbour B-B correlation is strongly antiferromagnetic. The corresponding correlations in 
the quantum model are both different from zero and are of the same order of magnitude. 
One could argue that quantum fluctuations favour planar versus non-planar arrangement of 
spins. This argument is supported by (i) the circumstance that there is a planar classical 
state of almost the same energy as the non-planar state having finite nearest-neighbour B-B 
correlations and (ii) by investigation of the so-called scalar chirality Wijk = Sj • (Sj x S&) 
being nonzero only in non-planar states. This kind of order parameter was widely discussed 
for the Ji — J 2 model [14,16]. We choose j G B as the running index and consider for i, j, k 
sites forming an equilateral triangle like sites 1, 23, 24 in Fig. 4, i.e. we have k = j + &2 and 
i = j + ai + |ai. Then we use as order parameter (cf. [16]) 



where tj is a staggered factor being +1 on sublattice B x (i.e. sites 17, 19, 21, 23 in Fig. 
4) and -1 on sublattice B 2 (i.e. sites 18, 20, 22, 24 in Fig. 4). As shown by Fig. 8 this 
chirality is indeed large in the classical nonplanar states but we do not see significantly 
enhanced chiral correlations in the quantum ground state. 

Further increasing J AB leads to a transition from the complex singlet S = phase directly 
to an S = 2 phase at Jab ~ 3.13, which is the quantum counterpart to the classical canted 
ferrimagnetic phase IV. This transition is very close to the classical transition III-IV which 
is also a transition from zero S to tai to finite S to tai- 

The last transition is that to the collinear ferrimagnetic state with S — 4 at Jab ~ 3.61. 
This value is significantly smaller than the corresponding classical value, again indicating 
that quantum fluctuations favour collinear spin ordering leading to an enlarged stability 
region of the collinear ferrimagnetic phase. Notice, that the additional jump just before the 
last transition is attributed to a change in total spin S from S = 2 to S = 3 corresponding 
to the increase of S to tai m the classical phase IV. One characteristics of both ferrimagnetic 
(S > 0) phases are the positive correlations within a subsystem (see Figs. 5 and 6) but 
negative correlations between the subsystems (see Fig. 7). 




(9) 
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IV. LINEAR SPIN- WAVE THEORY FOR THE NEEL PHASE 



The parameters for which the Neel phase I is realized most likely correspond to the situa- 
tion in (Ba, Sr) 2 Cu^OiCl2- Therefore we present a more detailed analysis of the magnetic 
ordering of this phase using a linear spin-wave theory. Within this approach we calculate 
the excitation spectrum, the order parameter as well as the spin-wave velocity. 

As usual we perform Holstein-Primakoff transformation. Because the magnetic unit cell 
contains six spins we need at least six different types of magnons being distinguished by a 
running index as illustrated in (Fig. 2). After transforming into the k-space the Hamiltonian 
(1) reads 

H = -ANJ AA s 2 - 2NJ BB s 2 + £ # k 

k 

with 



#k = ^Jaas [aj k a iyi + aj k a 2k + a3 k a 3k + aJ k a 4k J + AJ BB s (ag" k a 5k + ag" k a 6k J 

- Jaas (7ik72k + 7*k72k) ( a ik a 4-k + a 2k°3-k + a ik0 4 -k + a 2 ka 3 -k) 

- Jaas (7ik72k + 7*k72k) (aik a 3-k + a 2k a 4-k + a ik0 3 -k + a 2 ka 4 -k) 

- Jbbs (7i k + 72k + 7it + 72k) (4k a t~k + ask^e-k) 

+ Jabs (g + l) (71k (a| k a 6 k + a 2k a6 k ) + 7ik ( a i-k a tk + a^ek)) /2 

+ Jabs (g - 1) (71k («ik a 6-k + o 2 ka 6 -k) + 7ik ( a ika 6 -k + a 2k a 6-k)) /2 

+ Jabs (g + l) (72k (a| k a5k + a 3k a^ k ) + 72k ( a 4ka^ k + aj" k a 5k )) /2 

+ Jabs (g - 1) (72k (a4 k a^_ k + a 3 kfl5-k) + 72k (o4ka 5 -k + a 3k a^_ k )) /2 

- Jabs (g - 1) (71k (a2k°5k + ai k a^ k ) + 7* k (a 2 kaj k + af k a 5k)) /2 

- Jabs (g + 1) (71k (<4k a 5-k + a ika 5 -k) + 7ik ( a 2ka 5 -k + ai"k a 5-k)) / 2 

- Jabs (g - 1) (72k (a 3k a 6k + a 4k a^ k ) + 7^ (a 3 k<4 k + a+ k a 6k )) /2 

- Jabs (g + 1) (72k (a 3k c4l k + 04ka 6 -k) + 72k ( a 3ka 6 -k + a4k a e-k)) / 2 



(10) 



and 7„ k = exp (ika„/2). Here N is the number of geometrical unit cells N = Af/3. The 
vectors a n are the unit vectors of the geometrical lattice: a.± = a (1, 0) and a 2 = a (0, 1) and 
a is the lattice constant (see Fig. la), g is defined as g — cos (</?), where ip parametrizes 
the angle between the Neel states of the classical subsystems A and B. Without any further 
calculation it is obvious that quantum fluctuations stabilize collinear ordering. According 
to the Hellmann-Feynman theorem [19] the relation dE/dX = (dH/d\) holds, where if is a 
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Hamiltonian depending on a parameter A and E is an eigenvalue of H . Because (10) depends 
on cos (</?)-terms only one finds dE/df ~ sin (<p) being zero for (p = 0, n, i.e. as discussed 
already above in the quantum system the classical degeneracy is lifted and collinear spin 
structures are preferred. Thus, all quantities have to be calculated as averages over both 
possible ground states belonging to ip — and n. 

The diagonalization of the bosonic Hamiltonian is carried out as usual by means of Green 
functions. As it should be there are six non-degenerated spin-wave branches - two of them 
are optical whereas the remaining ones are two acoustical branches per subsystem. The 
acoustical branches become zero in the center of the Brillouin zone, only. Expanding these 
branches in the vicinity of k = gives two different spin- wave velocities ca and cb- 



ca = as^2J\ A + U BB - -^f^ + q , c B = as^2J\ A + U BB - _ q 



gj2 j4 

\ ( 2 ^1a ~~ 4J| B ) + —2 — (J aa Jbb - JaaJbb) + ~pr~ (Jbb ~ Jaa) > (H) 

\ J AA J AA 



where the two acoustical branches belonging to the same subsystem have identical spin- 
wave velocities. At the classical phase-transition point Jab = 2V JaaJbb cb becomes zero 
wheras ca remains finite. 
The ground-state energy E™ is given by 

6 

E™ = -2N (2 J AA + Jbb) s (s + 1) + ]T £ u mk /2 (12) 

k m=l 

and the sublattice magnetizations is calculated by 

2 

(Sn) = S ~ — («nk«nk) , U = 1, . . . , 6 (13) 
iV k 

for A-spins as well as for _B-spins. 

The results of the spin- wave calculation for the relevant parameters Jaa = 1 and J B b = 0.1 
are presented in Figs. 9 and 10. We start with the ground-state energy shown in Fig.9. 
While the classical energy in phase I is independent of Jab we find a slight decrease with 
Jab in the quantum model. For comparison we show the exact-diagonalization and the 
spin-wave results for Af = 24. The difference is small (1.3% for Jab = 0) indicating that 
linear spin-wave theory seems to be well appropriate for phase I. 

The spin-wave theory allows to calculate the corresponding sublattice magnetizations (S A ) 
and (S B ) in the A and B subsystems for Af — > oo (see eq. 13). The results are shown in 
Fig. 10. Although (S A ) is slightly diminished with growing Jab the Neel order of the A 
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subsystem is stable within the limits of the classical phase I. Contrary to that the Neel order 
of the B subsystem is stronger suppressed and breaks down at Jab ~ 0.58. This finding 
is supported by the ED results for the spin-spin correlations (see Fig. 11), where we see 
also a stronger suppressing of B-B correlations with growing Jab than of A-A correlations. 
Hence we argue that the strong quantum fluctuations in the spin-half model in combination 
with strong frustration may lead to a novel ground-state phase with Neel ordering in the 
A subsystem but quantum disorder in the B subsystem. A similar observation recently has 
been made for the frustrated square-lattice J\ — J2 spin-one spin-half ferrimagnet, where for 
strong frustration the spin-half subsystem might be disordered but the spin-one subsystem 
is ordered [11]. 

V. SUMMARY 

In this paper the results of exact diagonalization and linear spin-wave theory for the ground 
state of a system of two interpenetrating spin-half Heisenberg antiferromagnets on square 
lattices are presented. We consider intra-subsystem couplings of different strength J BB = 
O.IJaa which correponds to the situation in (Ba, Sr^Cu^OiCh- In addition to strong 
quantum fluctuations there is a competing inter-subsystem coupling Jab between both spin 
systems giving rise to interesting frustration effects. The classical version of our model 
possesses a rich magnetic phase diagram with collinear, planar and non-planar ground states. 
Quantum fluctuations may change the ground-state phases. In particular, we find indications 
for preferring collinear versus non-collinear and planar versus non-planar phases by quantum 
fluctuations. 

For small Jab both spin subsystems are in the Neel state. These Neel states decouple clas- 
sically. Quantum fluctuations lead to a fluctuational coupling of both subsystems. With 
increasing Jab the frustration tends to destroy the Neel ordering of the weaker coupled B 
subsystem but not in the stronger coupled A subsystem. The comparison between exact 
finite-size data and approximate spin-wave data gives a good agreement between both ap- 
proaches. 
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TABLES 



phase 


range of stability 


energy Jfa 


5a 
s2JV 


Sb 
sN 


Stotal 

sN 


I 

II 
III 

IV 

V 


< Jab < 2V JaaJbb 
2\/JaaJbb < Jab < ^V^Jaa 
2V2Jaa < Jab < (7) 
(7) < Jab < (8) 
(8) < Jab 


—\Jbb — \ Jaa 

1 J AB 4 t 

6 Jll ~ 3 J ^4 

1 4 T 

i^/y from eq. (6) 
— 3«7aB + §</rB + f </A4 






eq. (3) 
1 






eq. (3) 
1 






eq. (3) 
1/3 



TABLE I. The five classical ground-state phases for Jab,Jbb > and Jbb < Jaa, where 
Sa(b) = I J2neA(B) Si I is the total spin of subsystem A(B), and Stotal = I J2 n &i\ 1S the total spin of 
the whole system. 
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FIGURES 
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a. 




Jaa 
Jbb 



J 



AB 



a 



FIG. 1. a: The A-B lattice with its geometrical unit cell defined by the unit vectors a x = (a, 0) 
and &2 = (0, a). Empty (filled) circles denote A(B)-spins, respectively, 
b: Exchange couplings of the model (1). 




FIG. 2. The classical ground state for Jab < %V JaaJbb (Neel phase I, cf. table I). Both 
subsystems possess Neel order. Because of the vanishing classical mean field both subsystems 
decouple magnetically and can freely rotate with respect to each other. The angle <p is parametrizing 
this degree of freedom. The magnetic unit cell containing six spins is given by t>i(2) = ai ± a.2- 
The spins within the unit cell are labeled by a running index n = 1, . . . , 6 corresponding to the six 
different magnons to be introduced in spin-wave theory. 
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FIG. 3. The spin orientations of phases II, III and IV. In the non-planar phase II there is an 
out-of-plane z component of the A spins illustrated by open and dashed circles, where the open 
circles and the dashed circles belong to opposite directions of the z component. In the planar phase 
III this out-of-plane z component is zero. In the non-planar phase IV the out-of-plane z component 
of A spins as well as of the B spins is uniform but opposite to each other. The magnetic unit cell 
of the magnetic states of phases II, III, IV contains twelve spins. 




FIG. 4. The finite A-B lattice with M = 3iV = 24 spins (periodic boundary conditions). The 
large circles (sites 1, 16) belong to subsystem A and the small circles to (sites 17, 24) to 
subsystem B. 



13 




1 2 3 4 5 



J AB 

FIG. 5. A-A spin correlation (SfSf) for = (1,2) and = (1,7) (see Fig. 4) for the 
classical (thin lines, length of classical spin vectors is choosen as s = 1/2) and the quantum (thick 
lines) model (J\f = 24, Jaa = 1, Jbb = 0.1). 
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, , for = (17,24) and = (17,19) (see Fig. 4) for 



FIG. 6. B-B spin correlation (S? S B 
the classical (thin lines, length of classical spin vectors is choosen as s 
(thick lines) model (M = 24, J AA = 1, J B b = 0.1). 
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FIG. 7. A-B spin correlation {SfSf ) for (i, j) = (1, 17) and (i, j) = (1,24) (see Fig. 4) for the 
classical (thin lines, length of classical spin vectors is choosen as s = 1/2) and the quantum (thick 
lines) model (M = 24, J AA = 1, Jbb = 0.1). 




FIG. 8. Scalar chirality W (see eq. (9)) for the classical (solid line, length of classical spin 
vectors is choosen as s = 1/2) and the quantum (crosses) model (Af = 24, Jaa = 1, Jbb = 0.1). 
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FIG. 9. Ground-state energy per spin for s = 1/2: spin-wave (Af = oo and Af = 24) and 
exact-diagonalization results and (Af = 24). 
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FIG. 10. Sublattice magnetizations (S\) and (Sg) (see eq. (13)) for s = 1/2 in the infinite 
system: spin-wave results (Jaa = 1, Jbb = 0.1). 
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FIG. 11. Nearest-neighbour correlation within the A subsystem and the B subsystem (exact 
diagonalization results for N = 24, Jaa = 1> Jbb = 0.1). For better comparison we have scaled 
(SjSj) by its corresponding values for Jab = 0. 
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